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Abstract 

The geometric approach to autonomous classical mechanical systems in terms of a 
canonical first-order system on the Whitney sum of the tangent and cotangent bundle, 
developed by R. Skinner and R. Rusk, is extended to the time-dependent framework. 



1 Introduction 

In 1983, it was shown by R. Skinner and R. Rusk that the dynamics of an autonomous 
classical mechanical system, with configuration space Q, can be properly represented by a 
first-order system on the Whitney sum T*Q(BTQ [jl], |^, If the system under consideration 
admits a Lagrangian description, with Lagrangian L G C°°{TQ), the corresponding first- 
order system on T*Q®TQ is a Hamiltonian system with respect to a canonical presymplectic 
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structure. The Skinner-Rusk formulation can be briefly summarized as follows. Denoting 
the projections of T*Q ®TQ onto T*Q, resp. TQ, by pri, resp. pr2, and putting oj = prlug, 
with LjQ the canonical symplectic form on T*Q, one can consider the following equation 

izuj = dH, (1) 

where H := {pri,pr2) — L o pr2, and ( , ) denotes the natural pairing between the dual 
bundles T*Q and TQ. If the given Lagrangian L is regular, analysis of (|l]) shows that there 
exists a unique solution Z which is tangent to the graph of the Legendre map Legi : TQ — > 
T*Q, (g^,t>^) I— > {q^,dL/dv^), where the are local coordinates on Q and {q^,v^) denote 
the corresponding bundle coordinates on TQ. 

Remark. Here and in the sequel we will adopt the following definition for the graph of a 
bundle mapping. Let Ei and E2 denote two fibre bundles over the same base space M, and 
let f : El E2 , (m, e) 1— > (m, fm{e)) be a fibre bundle mapping over the identity. Then we 
define the graph of / as the image of the mapping f Xm idsi '■ Ei E2 x m Ei , (m, e) G 
El h-> (m, /^(e),e). 

The equations of motion induced by the vector field Z are equivalent to the Euler-Lagrange 
equations for the system under consideration (see 0). The important point now is that 
this equivalence between a Lagrangian system and the corresponding first-order system ([l|) 
on T*Q © TQ also holds if the given Lagrangian L is singular. In that case, in order to 
extract a consistent system of differential equations from (|l]), one will have to invoke a 
constraint algorithm. In fact, one of the main motivations for the work of Skinner and Rusk 
was precisely to show that the Dirac-Bergmann approach to singular Lagrangian systems 
can be properly described on the Whitney sum of the tangent and cotangent bundle of the 
configuration space, thereby avoiding some ambiguities occurring in the literature on the 
subject. 

The Skinner-Rusk formalism has been applied by several authors in various contexts |^, ^ 
. As one of the benefits of the formalism it turns out that it provides an appropriate setting 
for a geometric approach to constrained variational optimization problems. The latter are 
frequently encountered, for instance, in mathematical economics and in engineering. This 
has been demonstrated, in particular, for some optimal control problems in Q, and for the 
so-called vakonomic dynamics in |^ (see also [^). The fact that it would be interesting to 
extend those results to the time-dependent framework, allowing for systems with an explicit 
time-dependence of the 'forces' and/or the constraints, is the main motivation underlying the 
present work. More precisely, we will develop here a time-dependent version of the Skinner- 
Rusk formulation of dynamics, using the language of jet bundle theory [0] and cosymplectic 



geometry |jTO[ . Among the virtues of this new formulation of time-dependent mechanics, we 



would like to stress the possibility it offers to model a large class of systems, also in areas 
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such as economics and control theory. Apphcations for which this framework seems to be 
particularly well-suited include stabilization and trajectory tracking of mechanical systems 



by means of time-dependent transformations (see e.g. |TT 



Our starting point is a fibre bundle tt : E — > M, with E representing the evolution space of a 
mechanical system. Although it is quite common in treatments of time-dependent mechanics 
to fix a trivialization of vr, i.e. to work on a direct product space M x Q{= E), we will not 
select such a trivialization here. The natural space to consider then for the treatment of 
time-dependent Lagrangian mechanics is the first jet space J^vr, with the Lagrangian of the 
system being given as a function L G C°°(J^7r). The immediate candidate for replacing 
the direct sum T*Q © TQ in the Skinner-Rusk model for autonomous systems, seems to be 
the fibred product J^vr* x e J'^tt, where J^vr* is the 'dual' of the affine bundle J^tt (for this 
notion of dual, see e.g. It turns out, however, that this is not an appropriate choice 

for the following reasons. First, there does not exist a natural pairing between J^ir and 
J^TT*, needed for the construction of the time-dependent analogue of the "Hamiltonian" Ti. 
appearing in (|l]) and, secondly, there is no canonical 2-form on J^tt* to take over the role 
of the symplectic form ujq in the autonomous picture. To overcome these difficulties we will 
show that the appropriate space to consider is the fibred product T*E x e J^t^- 

Two final remarks are in order here. First of all, although we will restrict ourselves to 
Lagrangian systems, it is clear that, in analogy with the autonomous case (cf. the 
treatment can be easily extended to more general time-dependent mechanical systems, with 
forces not necessarily derivable from a potential. Secondly, it is interesting to note that the 
ideas developed here also admit a further extension to classical field theory, as has been 
demonstrated in a recent paper by M. de Leon et al. |T3I . 



The present paper is organized as follows. In the next section we briefiy recall the jet 
bundle approach to time- dependent Lagrangian and Hamiltonian mechanics, including the 
description of the constraint algorithm in case of a singular Lagrangian. In Section 3 we 
then develop the Skinner-Rusk formalism for time- dependent systems. We end with some 
conclusions and an outlook on future work along these lines. 



2 Non-autonomous Lagrangian systems 
2.1 The regular case 

Let TT : E — > M be a fibre bundle (the evolution space), with dimE = n + 1 and local 
bundle coordinates {t,q^),A = 1, . . . Consider the corresponding 1-jet space J^vr, with 
coordinates (t, q"^, q^) and associated projections tti : J^vr — > M and tti^q '■ J^t^ — ^ E. 
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Given a time-dependent Lagrangian L : J^n — > M, the Euler-Lagrange equations read in 
local coordinates 

These equations can be rewritten in geometrical terms as follows. Define the Poincare-Cartan 
1-form and 2-form 

eL = Lfj + s''{dL), ujL = -deL, 

d 

where fj = rrUdt) and S = idq^ — q^dt) ® 7— is the canonical vertical endomorphism on 

J^n (see 0). The action of S on 1-forms is denoted by S*. The equations (^ can then be 
expressed as 

ix^^L = , = 1 • (3) 

If the given Lagrangian is regular, then ujl has maximal rank and the pair {ujL,fj) determines 
a cosymplectic structure on J^vr, i.e. both forms are closed and satisfy the conditions u'^Afj 7^ 
0, ul^^ = (cf. ig). It then follows that (g) admits a unique solution, called the Euler- 
Lagrange vector field for L, and which we will denote by Xl- It is a second-order vector 
field, i.e. S{Xl) = and ixLV = 1, and a direct computation shows that integral curves of 
Xl determine solutions of the Euler-Lagrange equations (|^) and vice-versa. 

There also exists an alternative Hamiltonian description of the problem. Consider the Leg- 
endre map Lcql '■ J^t^ — ^ T*E, defined by LegL{j}(f)){v) = {QL)ji(t,{v) for j^cj) G J^n, 
V G T^{t)E and for any v G Tji^J^n such that TTifi^{v) = v. Let Vn denote the subbundle 
of TE consisting of vr-vertical tangent vectors, and denote its annihilator in T*E by (Vti)". 
Consider the quotient bundle J%* = T*E/{yTxY, which is called the dual of J^vr, with 
associated projections v : T*E — > J^tt*, tti q : J^n* — > E and tti : J^tt* — > M. Finally, 
denote by legi '■ J^t^ — ^ J^vr* the composition legi = u o Leg^. If L is regular, then LegL 
is an immersion and legi is a local diffeomorphism. 

Assume now that the Lagrangian L is hyperregular, that is, legL is a global diffeomorphism. 
Consider then the map h = Legi o leg^^. The mapping h : J^n* — > T*E is a section of the 
projection u, i.e. v oh = id^^ , and is called a Hamiltonian of the system. Next, denote by 
uoe the canonical symplectic two-form on T*E and let uj^ = h*ujE be its pull-back to J^tt* 
under h. If rji := {niYdt, then {LJh,Vi) defines a cosymplectic structure on J^tt*. In addition, 
one has that legKuh) = ojl and leg^rji = fj. It then easily follows that the solution X of (g 
is legL-related to the (unique) solution Y of the equations 

iyUJh = , iyVi = 1 • (4) 
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Note that, always under the assumption of (hyper-) regularity of L, the vector fields Xl, resp. 
F, are precisely the Reeb vector fields corresponding to the cosymplectic structures (0;^,//), 
resp. {uJh.rji). 



2.2 The singular case: the constraint algorithm 

Suppose now that the given Lagrangian L is degenerate, in the sense that the Hessian matrix 
{d^L/dq^dq^) is singular. We confine ourselves to the case where this matrix has constant 
rank everywhere, say r. The pair {uL^rj) then satisfies the following relations (cf. fl^, |15[): 

It follows that 2r < rank uol <2r + 2. In general, the equations (^) will not admit a global 
solution X. Moreover, if a solution exists it will not be unique. Therefore, in order to 
determine a consistent dynamics for such a system (if it exists) one has to apply a constraint 
algorithm which, at least in the case of Lagrangian mechanics, should be supplemented 
with the "second-order equation condition". For completeness, we will now briefly sketch 



the constraint algorithm described in |Tj, ^ , which is an adaptation to the time- dependent 
setting of the well-known geometric constraint algorithm for presymplectic systems developed 
by M. Gotay and J. Nester P, [E 



With a view on its application later on, we will describe the constraint algorithm here in 
the general framework of a structure (M, f2, rf) consisting of a smooth manifold M, a closed 
2-form Vt and a closed 1-form r], satisfying 



A 7/ ^ , ^r+' Ari = o , ^r+' = o , 

for some r < dimM. On M we then consider the equations 

ix^ = , ixV = '^ ■ (5) 

One can prove that there exists a vector X^. G T^-M satisfying these equations at the point 
X iff rank = 2r (see [|13)- In particular, it follows that (^) admits a global (but not 



necessarily unique) solution X iSQ has constant rank 2r, in which case the given pair {Q, rf) 
defines a so-called precosymplectic structure on M. If this is not the case, the constraint 
algorithm proceeds as follows. Put Pi := M and consider the set 

P2 := {x e M \ 3X^ e T^M such that ix^n^ = , ix^Vx = 1} • 

According to the previous observation, this set can be equivalently characterized by P2 = 
{x G Pi I rank Qr^ = 2r}. We then assume that P2 is an embedded submanifold of Pi{= M) 
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and we denote the natural inclusion by j2 : P2 ^ Pi- We are then assured that the equations 
(I) admit a solution X defined at all points of P2, but X need not be tangent to P2 and, 
hence, does not necessarily induce a dynamics on P2. We therefore have to continue the 
constraint algorithm by considering the subset 

P3 : = {x G P2 I G T^P2 such that ixM^) = , ixM^) = 1} 

= {xeP2\vxe KT.P2)} , 

where b is the bundle morphism defined by 

\?:TM T*M , V G T,M ^ + • 

Assuming P3 is a submanifold of P2, with inclusion map ja : P2 "—>■ Pi, it follows that there 
exists a vector field X on P2, which satisfies (j^) at points of P3. Again, however, such an X 
need not be tangent to P3, and one may have to repeat the above procedure. In this way, a 
descending sequence of submanifolds 



'<^^P,c^...^P3^P2^Pl(= M) 



is generated, with 



P, := {x G P,.i I 7], G b(T,P,_i} > 2) , (6) 

and where Pe is called the £-ar?/ constraint submanifold. If this sequence terminates at a 
nonempty set, in the sense that for some finite > 1 we have Pk+i = Pk, but P^ 7^ Pfc^i, 
then Pfc is called the final constraint submanifold, which we denote by Pf. Now, it may still 
happen that dimPj = (i.e. Pf is a discrete set), in which case the given problem admits no 
proper dynamics. However, if dimPj > 0, then we know by construction that there exists a 
vector field X on M, defined along Pf, which is tangent to P/ and satisfies the equations 

ix^Pf = , ixV\Pf = 1 5 

i.e. the given dynamical problem admits a consistent solution on Pf. In general, however, 
this solution will not be unique: given a solution X, for any smooth section Y of the bundle 
(ker Q fl ker rj) fl TPf over Pf, X + Y is also a solution. 

If we are dealing with a time-dependent Lagrangian system, i.e. with M = J^tt, Q = ujl, rj = 
fj {= dt), this is not the full story. First of all, we then also have to impose the so-called 
"second-order differential equation" condition, i.e. we are only interested in a solution X 
which determines a system of second-order ordinary differential equations. Secondly, as in 
the autonomous case, one can develop a similar constraint algorithm on the Hamiltonian 
side (i.e. on JV*) and, under a suitable assumption regarding the given Lagrangian, one 
can show that both descriptions are equivalent. For more details, we again refer to [14, 15]. 

In the next section we will show how the above can be translated into a Skinner-Rusk type 
formulation for time-dependent Lagrangian systems. 
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3 Skinner-Rusk formulation 



We start again from a fibre bundle n : E ^ with n- dimensional fibre, and its first jet space 
J^TT. Bundle coordinates on E and J^n are denoted by {t,q^) and {t,q^,q^), respectively. 
Canonical coordinates on T*E will be written as (t,q^,T,pA) and the canonical symplectic 
form on T*E then reads cue = dq"^ AdpA + dtAdr. We now consider the fibred product of the 
bundles T*E and J^n over E, i.e. T*E x e J^t!', with projections pri : T*E x e J^t^ T*E, 
pr2 : T*E x^; J^vr J^vr and pr : T*E Xe J^tt E. The natural bundle coordinates on 
T*E X E J^TT are {t, g^, r,p^, g^). 

On T*E X J^TT we define the 2-form uj as the pullback of the canonical symplectic form 
on T*E, i.e. u = pr\ujEi ^"^^ the 1-form rj = {tt o pr)*{dt) = pr^^fj. For simplicity we will 
sometimes write r] = dt. Recall that the affine bundle J^tt can be identified with an affine 
subbundle of TE whose underlying set is given by {v G TE \ {dt,v) = 1}. In coordinates, 
the natural embedding j : J^n ^ TE reads j(t, g^, g^) = {t, g^, 1, q^). 

Given a Lagrangian L G C°°(J^7r), we can define the following function on T*E Xe J^tt'- 

'H = {pri,j o pr2) - pr^L , 

where ( , ) denotes the natural pairing between vectors and covectors on E. In coordinates 
this becomes TC = pAq"^ + r — L(t, g^, g"^). Putting 

ujt-c = u + dlH. A Tj , 

we can then consider the following equations: 

izUJH = 0, izr] = '^- (7) 

Let us try to find out, in coordinates, what kind of dynamics is encoded by (|^). For that 
purpose, we write Z as 

d d d d d 



dt ' " dq^ "dr """dpA ' """dq' 
From the second equation in (|^) we deduce that Zt = I, and the first equation then becomes: 

iz^H = 'iz^ + ZiTi.)dt — dH 

f)r f)T r)T 

= {Z{n) + Zr)dt + - Z.Jdq^ - {pA - ^)dq^ + {Z^A - q^)dpA 

= 0. 

dL 

This immediately gives Z^a = g^ and Zp. = 7—7, together with the following constraint 
dL 

equations: pa = 'q^- These constraints determine a submanifold of T*E x e J^tx which, for 
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convenience, we will denote by M^. With the above expressions for Z^a, Zp^ and Zt, we see 

dL 

that the remaining condition ZiTi) + — = is identically satisfied at all points of Ml, 

ot 

irrespective of the value of the components Z^ and Z^a. Note that the relation ZqA = 
refiects the second-order differential equation property. 

The previous analysis already shows that (|^ locally admits a solution Z, defined in points 
of the submanifold Ml of T*E x e J^tt- In fact, we have a whole family of solutions since 
the components Zr and Z^a can still be chosen arbitrarily. 

In order to obtain consistent equations of motion, however, we have to impose the condition 

dL 

that Z be tangent to the submanifold Ml, that is, the functions Z(pA — tt^t) should vanish 
at points of Ml for all A = 1, . . . ,n. We now have that 

dL, dL d^L .R d^L ^ d^L 



QgA QgA dtdq^ dq^dq^ ^ dq^dq^ 

Clearly, if L is regular, the vanishing of @ fixes all the components Z^a as functions of 
(t, g^, g^) on Ml- If not, one will have to apply a constraint algorithm. 



3.1 The regular case 

Let us assume that the given Lagrangian L is regular. The previous analysis tells us that the 
system admits a solution Z on Ml and it follows from the expression for the components 
Zt, ZqA, Zp^, ZqA that any integral curve of Z on Ml determines a solution {q"^{t)) of the 
corresponding Euler-Lagrange equations of motion 

d f dL\ dL 



0, {A = l 



dt \dq^ J dq^ 

However, the solution Z is not unique since the component Z^- is still undetermined. This, of 
course, is not surprising since d/dr belongs to ker ujt-i fl ker i], i.e. ia_ujn = and id_ri = 0. 
In order to obtain a unique dynamics, we now impose an additional constraint 

T = L-q^ 



dq 



dL 



Together with the constraints pa = these are (locally) the defining equations of a 

submanifold of T*E x ^ J^tt, namely the graph of the (extended) Legendre map 

LegL : T*E , {t, q^, q^) ^ {t, L - """^ ^ 



dq^ ' dq 



(for the intrinsic definition of LegL, see e.g. [^, and for the notion of graph considered here 



see the Remark in the Introduction). We denote the graph of Lcql by graph^. Clearly, 
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graph C Ml and if we now require that Z should be tangent to graph^, it readily follows 
that 

Zr = ZiL - g^^) , (9) 

which uniquely fixes Z^-. Note that the differential equation corresponding to the r-component 
of Z represents the so-called "energy-balance" equation from time-dependent mechanics. 

The above construction was carried out on an arbitrary natural bundle chart of T*E x ^ J^tt. 
Using a standard argument it then follows that Z is in fact well-defined on the whole of 
graphj^. Defining the mapping 



LegL XEidji^ : J^tt — > T*E XeJ\, {t, q^, q^) i — ^ (t, g^, L - q' 



dL dL 



dq^ ' dq"^ ' 



we see that lm{LegL x e idj^-jr) = graph^;^ and it is not difficult to verify that the unique 
solution Z of (0), defined on graph^, and the Euler-Lagrange vector field on J^n are 
related by 

{LegL Xe idji^)^{XL)x = Z^ , 
where x = {LegL x e idjiT,){x), for all x G J^n. 

The previous discussion can now be summarized by the following proposition. 

Proposition 3.1 For a regular Lagrangian L, the system ^ admits a unique solution Z 
defined on, and tangent to, graphic, and the induced vector field on graphj^ is {LegL x Eidji^^)- 
related to the Euler-Lagrange vector field on J^tt. 

Recall that there exists a canonical projection u : T*E J^tt* {= T*E/{Vit)^) which, in 
coordinates, reads ^{t, g^, t,pa) = {t, q^,PA) (cf- Section 2.1 and [Q)- Let us assume that L 
is hyperregular such that the mapping legL '■= voLegL is a global diffeomorphism. Consider 
the fibred product J^vr* x e J^t^, with associated projections Ai and A2 onto J%* and J^vr, 
respectively. We can then define the following projection 

u XEidji^ : T*E Xe — > J^n* Xe J^n , {t, q^,T,pA, q^) 1 — > {t,q^,PA,q^) ■ 

From the discussion above we deduce that, along graph^, the vertical distribution determined 
by the projection u Xe idji^^, i.e. ker(z/ x e ^c^jItt)*, is invariant under Z in the sense that 

[Z, ker(z/ x^; idji^)^] C ker(z/ x^; idji^)^ . 
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Hence, in the hypcrregular case, the solution vector field Z is projectable onto J^tt* x e J^t^- 
Its projection locally reads 

This is the unique vector field on J^n* x e J^t^ determined by the equations 

where we recall that h — Lcql o leg^^ and ujh — h*ujE (cf. Section 2.1). 

In the general case, however, it is not possible to represent the dynamics of the non- 
autonomous problem corresponding to L by a first-order system on J^tt* x e J^t^- 

3.2 The singular case 

Returning to the beginning of this section, let us now assume that L is not regular. Observe 
that, with Lu := prlcuE, we have 

uj^{— uJu A uJu) — OJ /\0J -\- luj /\ dJi A 77 , 

and, in general, 

u;^ = u;'' + kLj'"^ AdHAv, VA; . 
Herewith, it is straightforward to check that the pair {cu-h, rf) satisfies the following relations: 

o;;^ A 77 7^0, u;;^+^A77 = 0, u;;^+^ = 0. 

Indeed, we have 

o;;^ A 77 = a;" A 77 = (-l)^^'^7i! o?g^ A • • • A dg" A dpi A • • • A dp„ A 77 7^ , 
uj^^ A 77 = (-1)^"^'2"^'^ {n + ly. dt A dq^ A . . . A dq"" A dr A dpi A . . . A dpn Ar] ^ , 
oj^"^ = u;"+2 + (n + 2) AdHAdt^O. 

This implies, in particular, that 27i < ranko;-^ < 2(71 -|-1), where we recall that dim£' = 7i-|-l. 

Putting Ml :— T*E x e J^t^, we can now apply the constraint algorithm described in Section 
2.2 to the triplet (Mi, o;-^, 77). First of all, we consider the set 

M2 = {x e Ml I 3 e 7;Mi such that iz,0JH{x) = , izM^) = 1 } , 
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which can be equivalently characterized by M2 = {x G Mi \iai\kujT-i{x) = 2n}. In local 
coordinates we find 

u;+\x) = u''+\x) + [n + 1) u'' A dn A 7]{x) 

= u;"+^(x) - (n + 1) cj" A "^c^^ A r]{x) + {n + A c/g^ A r]{x) 

= {-l)^^^'^^{n + ly.-^^dq^ A ■ ■ ■ A dq'" A dpi A ■ ■ ■ A dp^ A dq^ A ■r]{x) , 

such that X G M2 if and only if 



on , dL 



{Pa - T77j)\x = , A=l,...n. 



dq^ w dq^ 



Observe that M2 coincides with the submanifold introduced at the beginning of this 
section. By construction, there exists a vector field Z on Mi, defined along M2, which 
verifies equations (0) at points of M2. But in general Z will not be tangent to M2 and so we 
then have to proceed with the constraint algorithm by considering the set 

M3 = {x e M2 I 3 e T,M2 such that iz^oo-uix) = , iz^vi^) = 1} • 

Assuming that M3 is a smooth submanifold, there will be a vector field Z defined along M3 
and tangent to M2, satisfying (J^) at each point of M3. Continuing this way, we obtain a 
descending sequence of submanifolds of Mi that, in the favorable case, will stop at a final 
constraint submanifold Mf on which there exists a consistent solution of the given dynamical 
problem (cf. Section 2.2). The constraint submanifolds Mi can still be characterized in an 
algebraic way similar to (^), with the map b : TMi T*Mi being induced here by the pair 

As in the autonomous case, we thus see that the constraint algorithm for time- dependent 
singular Lagrangian systems can be properly developed in terms of the structure {T*E x e 
J^tTjU-HjT]). To complete the picture, we have the following result which shows that this 
description is equivalent to the standard one based on the structure {J^n, 00^,7]). 



Proposition 3.2 Let {Pe}e>i, resp. {M£}e>i, denote the sequence of constraint submani- 
folds generated by applying the constraint algorithm to {J^n, ul, fj), resp. (T*Ex eJ^'^i v)- 
Then, for each i = 1,2, ... , we have that ipi+i = pr2\Mi+i '■ ^ Pi is a surjective sub- 

mersion such that the following diagram commutes 
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Ml := T*E Xe J^tt 



32 



Mo 



Ms 




^3 



Pi := J^n 

32 

Po 



(where ji : Pe Pe-i and j[ : Mi M^_i are the natural embeddings of the respective 
constraint submanifolds). Moreover, if there is a final constraint suhmanifold Mj := Mk C 
Ml (for some k > 2) on which there exists a consistent solution Z of then Z projects 
under ip^ onto a solution of (Qj on the final constraint suhmanifold Pf := Pk^i in J^tt and, 
conversely, any solution of (Qj, defined on Pf, is the projection of a vector field on Mj which 
is a solution of (Qj. 

The proof of this proposition essentially relies on the following two facts. First of all, the 
projection pr2 '■ T*E x e J^t^ J^tt has the appropriate "almost regularity" properties, that 
is: (i) pr2 is a surjective submersion and (ii) the fibres of this submersion are connected 



submanifolds of T*E x e J^tt, being diffeomorphic to 
computation shows that 



dq^ A d 
dq^ A d 



PL 
dq^ 
PL 
dq^ 



+ dEi A dt 
dL 



. And, secondly, a straightforward 



dq^ 



A dt 



dL 
dq^ 



dq^ A dt 



2| can be easily 



and, clearly, we also have pr^fj = rj. Herewith, the proof of Proposition 
completed, following the same reasoning as in the autonomous case 0]. 

The solution generated by the constraint algorithm (if it exists) is not unique. On the other 
hand, we may observe that if the given dynamical problem (|^) admits a consistent solution 
Z on a final constraint submanifold Mf then, by construction, its projection onto J^tt will 
automatically verify the second-order equation condition along a submanifold of Pf. This is 
again in full analogy with the situation encountered in the autonomous case. 
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Next, assume that the given Lagrangian L G C°°(J^7r) is almost regular in the following 
sense: (i) Legi^J^n) is a submanifold of T*E, (ii) Legi, regarded as a map from onto 
its image, is a submersion with connected fibres, (iii) legj}{legL{x)) is a connected set for 
all X G J^TT. In |TH] it has been shown that, with these assumptions, one can develop 
a constraint algorithm on J^tt* which is equivalent to the one on J^vr. Again as in the 
autonomous case (see 0), one can demonstrate that a solution of the constrained analysis 
on J^TT* can be related to a solution Z of (|^), defined on the final constraint submanifold 
Mf. This connection is established by choosing a suitable (local) section a of the projection 
uopn : T*E X E J^T^ —>■ J^TT* and restricting Z to Im(o") fl Mj (recall that u is the canonical 
projection of T*E onto J^tt*). 



4 Conclusions 

We have developed a non- autonomous version of the Skinner-Rusk approach to (Lagrangian) 
mechanics and have shown that, both in the regular and in the singular case, this yields a 
first-order system on the fibred product T*E x e </^vr which encodes all the information of the 
dynamics of the system under consideration. This approach to time-dependent mechanics 
possesses the same virtues as in the autonomous case, such as the fact that the 'Hamiltonian' 
TC is defined without having to solve the relations pa = dL/dq^ for (some of) the velocities. 

Within the above framework for the description of time-dependent mechanics, there are 
several lines of investigation that seem to be worth pursuing, such as: the role and the 
nature of gauge transformations in the case of singular Lagrangians and the general study 
of symmetries of (time-dependent) mechanical systems. In addition, it would certainly be 
of interest to use this formalism for establishing a geometric formulation of optimal control 
problems with an explicit time-dependence, such as in the case of time- dependent vakonomic 
dynamics, thereby generalizing the work presented in 
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